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Abstract. The connection of Dark Matter (DM) to our particle physics model is still one of the 
open cosmological questions. In these proceedings I will argue that axinos can be successful Cold 
Dark Matter candidates in models with Supersymmetry and the Peccei-Quinn solution of the strong 
CP problem. If they are the Lightest Supersymmetric Particle (LSP), they can be produced in 
the right abundance by thermal scatterings and out of equilibrium decays of the Next-to-Lightest 
Supersymmetric Particles (NLSPs). Moreover if the NLSPs are charged, their decay could help us 
understand which is the LSP, e.g. between axino and gravitino. 

INTRODUCTION 

The nature of Dark Matter is still unknown today, but from simulations of structure for- 
mation we know it must be a neutral, Cold, collision-less (i.e. quite weakly interacting) 
and very long lived particle yj]. Unfortunately such a particle does not exist in the Stan- 
dard Model (SM): neutrinos are neutral and massive, but so light that they are at most 
Hot Dark Matter. Therefore we are obliged to look for DM candidates in models beyond 
the SM. One of the best motivated extensions, supersymmetry with R-parity conserva- 
tion, naturally gives us a stable massive particle, the LSP. To be DM such a particle has 
to be neutral and very weakly interacting, so usually only the neutralino or the grav- 
itino are possible LSPs. But if we invoke the Peccei-Quinn solution to the strong CP 
problem, a new multiplet has to be introduced, the axion multiplet. As long as super- 
symmetry is unbroken, this whole multiplet remains light, so that no supersymmetric 
mass parameter is allowed for it (contrary than for the Higgses). After supersymmetry 
breaking the fermionic component, the axino, obtains a mass, but it still could be the LSP 
and make a very good DM component. We will present in this talk a summary of axino 
CDM J3, S 0, El 0] and explore in particular if such particles can be produced in suffi- 
cient numbers to make up most of the DM and what that implies for the supersymmetry 
breaking parameters and collider searches. 



PRODUCING AXINOS IN THE EARLY UNIVERSE 

We briefly review here the two main mechanisms that produce axinos in the Early 
Universe. We concentrate here on the hadronic type of axion models, where is it expected 
that the axion multiplet does not interact directly with the SM multiplets and therefore 
the axino does not mix substantially with the standard neutralinos. In the other type of 
models, this mixing can be larger and the production is therefore enhanced. 



Thermal scatterings 



Any particle, even very weakly coupled, is produced in the thermal plasma by scatter- 
ings of the particles that are in thermal equilibrium. Axinos couple directly to the gluons 
and gluinos due to the axion "anomaly" coupling 

where <t> a is the axion multiplet containing the axino a, W the gluon vector multiplet 
containing the gluino g and the gluon G^ v and f a is the Peccei-Quinn scale of the order 

of 10 11 GeV due to axion physics |2D. So many scattering of the primordial plasma 
involving colored particles can produce axinos l . The axino number density is given by 
solving a Boltzmann equation of the type 

— +3Hn a = 52(<7(«' + J -^a + ...)v re i)ninj + 2_,(r(i^a + ...))ni (2) 

where we are neglecting back-reactions, that are suppressed by ria <C «/. 

At high temperature the 2-body scatterings dominate the production, since they con- 
tain a vertex given by the dimension 5 operator in eq. CO) and they show a characteristic 
linear dependence on T . So most of the axinos are produced at the highest temperature, 
the reheat temperature Tr, and the axino number density is proportional to Tr. Some 
of the two body scatterings are IR divergent due to the massless gluon propagator; in 
the thermal bath such a divergence is screened by the presence of a thermal gluon mass 
~ gT. In our computation we introduced such IR cut-off by hand [3]. A self-consistent 
procedure is instead to perform a full resummation of the Hard Thermal Loops as in ||9|]. 

At lower temperatures the decay terms start dominating and the number density is no 
more proportional to the reheat temperature, it depends instead on the supersymmetric 
parameters, in particular the gluino and squark masses ||4|]. 

Using the expression for the present axino energy density as 

w-a / m \ =0-72eV — — — , (3) 
s(T) V 0.2 J ' 

where s(T) = 2.89 x 10 3 (jim) 

cm 3 is the present entropy density, we can then obtain 
a bound on the reheat temperature as shown in Figure 1 . 



Out of equilibrium decays 

An axino population is also generated by the NLSP decay after it freezes out from 
the thermal bath. The heavier superpartners cascade-decay quickly into the NLSP (or 



The same happens also in the case of the gravitino, but with different vertex structure and scale 1 8]. 
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FIGURE 1. Maximal reheat temperature as a function of the axino mass obtained by requiring that the 
axino energy density is below the present DM density [3]. The difference between solid and dashed line 
is due to the inclusion of the decay term in the Boltzmann equation (0. In the gray area we expect the 
non-thermal production via out of equilibrium decays to be also substantial. 



very rarely in the LSP as we discussed in the previous section) while still in equilibrium, 
but the NLSP instead has a lifetime longer than the freeze-out time: in fact all the axino 
couplings are suppressed by the Peccei-Quinn scale f a ~ 10 11 GeV and so the NLSP 
lifetime is of the order of seconds or longer. Then the freeze-out process is unaffected 
and the decay takes place only much later as shown in Figure 2. 

In this case, the axino number density can be directly computed from the NLSP 
would-be-relic number density as 

AT = — CInisp. (4) 

mNLSP 

If the mass ratio is not too small, we still have a connection with the classical WIMP 
mechanism in case the NLSP is a neutralino or a stau. 

On the other hand, a couple of problems can arise if the decay happens too late: 

• Big Bang Nucleosynthesis can be spoiled by the energetic "active" particles pro- 
duced in the decay along with the axino: the strong limits on the injection of en- 
ergetic particles depend on the electromagnetic/hadronic nature of the produced 
showers and the decay time ifioll . In general such limits are weak for the ax- 
ino case since the NLSP lifetime (excluding a strong mass degeneracy) is below 
I0 2 s, but they can affect the region of small mass for both the neutralino and stau 
NLSP HH. 




FIGURE 2. Freeze-out of the NLSP and subsequent decay into axino. Due to R-parity conservation 
the number of axino produced in the decay is the same as the NLSP number. 



• Are axino from the decay cold enough to be CDM ? They are relativistic at 
production even if the NLSP is not and have a non-thermal spectrum: 

ma 2m 5 \g*{T dec ) J T dec 

where T dec is the temperature at the decay time. So the question is if they have suffi- 
cient time to cool down before structure formation begins. In [11] such constraints 
have been studied and the conclusions is that an axino mass at least of order of 
1 GeV is probably needed. 



AXINOS AND THE CMSSM 

Depending on the parameters and Tr, either production mechanism can dominate and 
produce sufficient axinos to explain the present DM density. In general either Tr is 
bounded as in Fig. 1 or the axino is so light to be a subdominant (warm or hot) DM 
component. In the last case in our scenario the axion Q] could be the DM. 

Assuming that the axinos are CDM and that the supersymmetric partners of the 
SM particles can be described by the Constrained MSSM, we can see which is the 
preferred parameter region depending on the production mechanism. In the CMSSM the 
superparticle spectrum and couplings are simply function of the SM and additional four 
parameters: the ratio of the Higgs v.e.v. tan/3, the gaugino and scalar masses mi/ 2 i m o 




FIGURE 3. Allowed parameter space for the case of dominant thermal production 1 5] . We have chosen 
here Tr = 200GeV, m c ~, = 1 GeV and f a = 10 GeV. The dark gray strip gives axino in the right abundance 
to explain all DM, while the lighter gray areas are excluded by LEP constraints or too large axino number 
density. The white area has too low axino density to explain DM. 



and the trilinear coupling Aq, which are universal at the GUT scale. The modulus of the 
/l parameter is fixed by radiative electroweak symmetry breaking and we will always 
consider the positive sign in the following. 



Mostly thermal production 

In the case of high Trh all the particles in the thermal bath can be treated as massless 
and so there is practically no dependence on the supersymmetry breaking parameters. 
On the other hand if we require the axino mass to be above 1 GeV, the reheat temperature 
has to be sufficiently low and comparable to the superparter masses, so the decay term in 
the Boltzmann equation become important and a strong dependence on the gluino mass 
appears also due to the squark-quark- axino coupling [4]. We have then that the allowed 
region is a narrow band in the gaugino mass parameter with a much smaller dependence 
on niQ as shown in Figure 3. Note that in this case the small gaugino mass region is 
excluded because there too many axinos are produced. The exact position of the allowed 
band though strongly depends on the chosen reheat temperature and it moves at larger 
m 1/2 for larger 7>. 
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FIGURE 4. Allowed parameter space for the case of dominant production via out of equilibrium NLSP 
decay |5]. We have chosen here Tr = 200GeV, m ( ~, = 1 GeV and f a — 10 GeV. The dark gray region 
gives axino in the right abundance to explain all DM, while the lighter gray areas are excluded by LEP 
constraints or too large axino number density. The white area has too low axino density to explain DM. 

Mostly out of equilibrium NLSP decay 

In the CMSSM the NLSP can be either the neutralino or the stau. The latter happens 
in the wedge with low mo, that is usually considered excluded if the stau is the LSP. In 
our case though the LSP is the axino and even the stau wedge is viable; in particular 
there is a wide region where the stau decay can produce the right amount of axino DM 
as we see in Fig. 4. On the other hand there is also a tiny strip in the neutralino NLSP 
region, analogous to the neutralino DM region, but shifted slightly due to the rescaling 
in eq. ©. Both regions are practically unaffected by BBN constraints contrary to what 
happens for the gravitino case ifl^l . 

HOW TO DISTINGUISH THE LSP ? 

If the axino is the LSP, very different signals could be found at colliders depending on 
the nature of the NLSP If the neutralino is the NLSP, the only way to find out that it is 
not DM is if the mass and cross sections turn out to give a too large neutralino number 
density or to be excluded by direct DM searches. Then we would have good reasons to 
say that the neutralino must be unstable, but to study its decay will be very difficult. 

If the stau (or another charged sparticle) is the NLSP instead, we will have the striking 
signal of an apparently stable charged heavy particle in the detector. In that case it will 



be clear that the LSP must be a very weakly interacting particle, but to know which one, 
we will need to measure and study the NLSP decay. In particular to distinguish between 
axino and gravitino, that can give similar NLSP lifetimes, we will need to measure the 
branching ratio and the angular dependence of the radiative decay in order to reach a 



Axinos with masses in the MeV-GeV range are good CDM candidates for low re- 
heat temperature: they can be produced either from thermal processes or from NLSP 
decay with the right abundance. Such scenario is analogous to the gravitino LSP one, 
but an axino LSP evades more easily BBN bounds, since the NLSP lifetime is shorter 
than 10 2 s. Therefore both neutralino and stau NLSPs are allowed in our case. 

Compared to neutralino DM, different regions of the CMSSM parameter space be- 
come allowed and preferred, in particular even heavier sparticles masses and a charged 
NLSP like the stau. An apparently stable charged particle would surely give a striking 
signal at the LHC and ILC and indicate that the neutralino is not the DM. On the other 
hand, disentangling between LSP candidates will require to stop such NLSPs and mea- 
sure their decay; in particular axino and gravitino will be distinguished if a sufficient 
number of radiative decay can be observed. 
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